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We investigate the dual superconductivity mechanism in the exceptional group G2. This is a cen- 
terless group (no 't Hooft flux vortices are allowed) and we check for the presence of a magnetic 
monopole condensate in the confined phase by measuring on the lattice a disorder parameter re- 
lated to the vacuum expectation value of an operator carrying magnetic charge. The behaviour 
of the disorder parameter is consistent with the dual superconductor picture. A first step of an 
analysis on the thermodynamical properties of the theory is conducted by mean of this operator. 



XXIVth International Symposium on Lattice Field Theory 
July 23-28, 2006 
Tucson, Arizona, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Dual Superconductivity in G2 group 



Guido Cossu 



1. Introduction and Motivation 

Understanding confinement is one of tlie main issues in non perturbative QCD. In this respect 
identifying the degrees of freedom responsible of quark confinement is essential. Many choices 
are investigated in literature, between them center vortices and monopoles (both proposed by 't 
Hooft) are the most popular. The deconfinement transition in pure SU(3) gauge theory is triggered 
by the spontaneous breaking of center symmetry, this could suggest the direct relation between 
confinement and presence of a non trivial center. It's therefore an interesting question whether 
centerless groups have a confinement-deconfinement phae transition. In fact the S0{3) = SU (2) /Z2 
has been extensively studied on the lattice (see for example [|l|, |2l, ^). Moreover the exceptional 
group G2 is without 't Hooft center vortices and is the simplest one, that's why we focused on it. 

Another point of view concerning the confinement mechanism is the so called dual super- 
conductor picture, the subject of this work. As in an ordinary superconductor magnetic fields are 
squeezed in flux tubes (Abrikosov filaments) connecting magnetic charges due to condensation of 
Cooper pairs, quarks could be kept together by strings of (chromo)electric flux tubes due to con- 
densation of magnetic objects. According to this picture, the QCD vacuum, is a condensate of 
magnetically charged fields confining (chromo)electrically charged particles. 

The presence of strings connecting static quarks has been proved in several papers. A way of 
probing the vacuum state looking for the Abelian magnetic fields has been proposed and studied by 
the Pisa group in a series of papers [^, |5|, ^, ^J. Abelian monopoles are exposed after a gauge fixing 
and a "monopole probe" is inserted in the lattice: the expectation value of such an observable, 
usually called jJ., should be exactly zero if magnetic symmetry is respected, a value different from 
zero is instead a clear signature of the presence of an Abelian monopole condensate in the vacuum 
of the theory. In practice simulations are done by inserting a monopole field in a time slice and 
measuring a related quantity, the so called p parameter defined by 

p^-^ln{^) = {S)s-{SM)sM (1-1) 

in which Sm is the action with the "monopole probe" insertion. So p is the difference of the action 
minus the monopole action weighted with the monopole action in the path integral. For details we 
refer to the cited articles. 

The aim of this work is to give a criterion of confinement in G2 group by mean of the disorder 
parameter jj.. Strictly speaking there are no strings in G2, a situation similar to full QCD where 
string breaking occours. Anyway in full QCD with A/y^ = 2 [0] the operator /i is nonetheless a good 
order parameter. So, in the working hypothesis that this is a good observable for confinement, we 
look for dual superconductivity in the exceptional group G2. 

1.1 G2 group 

We are now going to state some basic facts about the Lie Group G2. Being defined in mathe- 
matics as the group of automorphisms of the octonions it can be naturally constructed as a subgroup 
of the real group S0{7) which has 21 generators and rank 3. To the usual properties of SO(J) ma- 
trices 

detn = 1 D.^^ = (1.2) 



2 



Dual Superconductivity in G2 group 



Guido Cossu 



we have in addition another constraint 

Tahc = Tilef^da^eb^fc (1-3) 

where Tabc is a totally antisymmetric tensor whose nonzero elements are (using the octonion basis 
given by g]) 

^123 = 7^176 = ^145 = T251 = 7246 = ^347 = 7365 = 1- (1-4) 



Equations 1.3 are 7 relations reducing the numbers of generators to 14. This is clearly seen also in 



this form which is actually used as a check during simulations: 

v = Q.v Vk = ^jiTijk k=\...l (1.5) 

The fundamental representation of G2 is 7 dimensional and by using the algebra representation 
of we can clearly identify an SU (3) subgroup and several SU (2) subgroups, with 6 of them 
we can cover the whole group, an important property for MC simulations. The first three SU (2) 
subgroups are the 4x4 real representation of the group while the remaining three are extremely 
difficult to simulate with standard techniques. See next section for details on simulations. 

The Lie group G2 has rank 2 as SU (3), this imphes that its Cartan subgroup, that is the maximal 
residual abelian subgroup after an abelian gauge fixing, is?7(l)x?7(l). Stable monopole solutions 
are classified according to the homotopy group' : 

K2{G2/U{\f) = Ki{U{\) X U{1)) = Z X Z (1.6) 

i.e. we have two distinct species of monopoles as for SU (3). 

Another interesting homotopy group shows that center vortices are absent in the theory: 

7ri(G2/^(G2)) = 7ri(G2)=0 (1.7) 

while for SU (3) for example 

;ri(5f/(3)/Z3)=Z3. (1.8) 

So G2 is a good playground to study the dual superconductor picture in a theory without center 
vortices, thus isolating monopole contribution in confinement. 



2. Simulation and Results 

Much work has been devoted to simulations on G2 Yang-Mills theory by the Bern group 
[ ]T0| , pT] ]. This group used the so called Fredenhagen-Marcu order parameter jr2| ] to demonstrate 
confinement in G2 theory, at least in the strong coupling limit. This parameter by construction 
makes sense only at zero temperature. Our proposal concerns a different approach using the mag- 
netic operator ji of use also in finite temperature simulations. Here we are going to investigate the 
thermodynamical and dual superconducting properties of the system. 

First some informations about the MC runs settings. To simulate the gauge theory 

^ = ;^TrF^v7>v (2.1) 
'First equality follows from 7t\ (G2) = 0. See for example PI 
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Figure 1: Plaquette susceptibility plotted against j3. The peak signals the bulk transition and the points 
enclosed by the square are the physical transition for the Ni — 6 lattices. The transition for A^, = 4 is too 
close to the bulk one and extremely difficult to measure (much of the signal is only "noise"). Also shown a 
simulation a T=0 on a 16'* lattice (blue points). 



we used a simple Cabibbo-Marinari update (heat-bath + overrelaxation in a tunable ratio, for every 
step) for the first three SU (2) subgroups (4x4 representation) spanning the SU (3) C G2. This 
simple setting cannot be used for the remaining three subgroups so we make a completely random 
gauge transformation every n updates (tipically 2) to guarantee the covering of the whole gauge 
group. Observables measured are the standard plaquette and the Polyakov Loop. To study the ther- 
modynamical properties we simulated several lattices of spatial dimension Ng = 12, 16,20,24,32 
and Nt = 6 (Nt = 4 only for the smallest lattice). At the transition we needed histories of more 
than 10"* updates. The code is very fast (using only real algebra), is written using explicitly SSE2 
instructions in single precision for the matrix-multiplication core and run on an Opteron farm here 
in the computer facilities of the Physics Department of the University of Pisa. 



2.1 Thermodynamics 

To guess where the physical transition should be, we looked at the plaquette and Polyakov 
loop susceptibilities. 

The first thing to notice is the presence of a bulk transition in the plaquette susceptibility 
(related to the specific heat), see fig. |TJ The peak is present in every volume and Nt always at 
the same j3 ~ 1.36. There's no scaling with volume and no movement toward the weak coupling 
region passing from A/^r = 4 to A/^t = 6 as we would expect for a physical transition. This big peak 
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Figure 2: Polyakov loop density (on the left) varying j3 in the range from 1.35, the critical coupling of the 
bulk transition "j3c", to 1.405 where we are in the deconfined phase (all data from the 6 x 16^ lattice). On 
the right a rescaling of the Polyakov loop for the smallest lattices assuming a first order transition (v = 1/3 
and yj^ = 2.47 based on a fit on the peak heights). The smallest lattice is evidently not in a scaling region 
while the other two agree quite well with the hypothesis. Points on the 6 x 16^ lattice on the left of j3c are 
probably bad sampled, we just need more statistics. 



completely overshadows the real physical transition that can be seen as a little peak in the weak 
coupling region at j8 ~ 1.395 for Nt = 6 and several spatial volumes. 

Polyakov Loop is insensitive to the bulk transition so we used it to detect the position of the 
physical one, even if it's not an order parameter, see fig. ^. From the behaviour of this observable 
one can guess a first order transition. An accurate scaling analysis is needed for a conclusive 
statement. 

2.2 Parameter p 

The p parameter which, we recall, is defined as the difference between the actions without 
and with the monopole insertion is expected to be volume independent as j3 — > and to show a 
dip at the physical transition, see for example [^] . Measures of this parameter showed it to be 
extremely sensitive to the bulk transition with a strong dip at jS ~ 1.36 for all volumes, see fig. ||. 
Dips at the physical transition are evident only for large volumes A^^ > 20, fig. ^. The background 
bulk transition precludes any direct scaling analysis using the monopole operator, moreover volume 
indepencence before the physical transition can't be seen; key features of the operator are obscured. 
Another method is needed to study the problem. 

A promising method for the p susceptibility measure is to renormalize the operator at T=0 
subtracting its value (po) = (p(r = 0)). This is analogous to a renormalization in the free energy 
density typical in thermodynamical analisies. This way is at an explorative stage. 

2.3 Conclusion 

We are working on the G2 centerless theory in search for confinement by dual superconduc- 
tivity. Our parameter p shows as expected a dip at the transition points but there's an unphysical 
background that obscures its scaling properties. We are developing a method to deal consistently 
with this background and find a first order scaling as expected from plaquette and Polyakov Loop 
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Figure 3: Results on the p parameter at different lattice volumes, strong dip occours at the bulk transition 
point. 
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Figure 4: Same as before, magnification of the physical transition region. 
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Monte Carlo hystories. Also a scaling analysis of the Polyakov Loop susceptibility, even if it is not 
an order parameter, seems to indicate a first order phase transition. 

We cannot yet conclude that the theory confines but first results, although plagued by the bulk 
transition background, seems to indicate a positive answer. 
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